Prepared for submission to JHEP 



Chiral perturbation theory for twisted mass QCD at 
small quark mass 



Satoru Ueda a and Sinya Aoki bc 

a KEK Theory Center, High Energy Accelerator Research Organization (KEK), 

Tsukuba 305-0801, Japan 
b Graduate School of Pure and Applied Sciences, University of Tsukuba, 

Tsukuba, Ibaraki 305-8571, Japan 
c Center for Computational Sciences, University of Tsukuba, 

Tsukuba, Ibaraki 305-8511, Japan 

E-mail: sueda@post.kek.jp, saoki@het.ph.tsukuba.ac.jp 

Abstract: We study the lattice cutoff (a) and quark mass dependences of pion masses 
and decay constants in the Nf = 2 twisted mass QCD, using the Wilson chiral perturbation 
theory to the next leading order (NLO). In order to investigate the region near zero quark 
mass, we introduce the power counting scheme where 0(a 2 ,am) terms are included in the 
tree level effective Lagrangian. At the NLO of this power counting scheme, we calculate 
the charged pion mass and decay constant as a function of the lattice cutoff as well as the 
twisted quark mass at the maximal twist. In this paper, we adopt two different definitions 
for the maximal twist. We confirm that the difference between the two appears as the 
0(a 2 ) effects so that the automatic 0(a) improvement is realized for both definitions. 
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1 introduction 

The twisted mass lattice QCD (tmlQCD) [1-3] has several advantages for numerical simu- 
lations, one of which is the automatic 0(a) improvement at the maximal twist [3-5]. The 
tmlQCD becomes free from 0(a) lattice spacing errors, by simply setting the twist angle 
to its maximum value u = w/2. 

This automatic O(a) improvement of the tmlQCD has been investigated in quench 
simulations by the XLF Collaboration [6-9] and Abedel-Rehime et al.[10, 11], while the 
unexpected first order phase transition has been found in full QCD simulations by DESY 
group [12-14]. Recently the European Twisted mass Collaboration (ETMC) has started 
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the large-scale full QCD simulations at the maximal twist with Nf = 2 [15, 16] and Nf = 
2 + 1 + 1 [17, 18]. 

The chiral perturbation theory (ChPT) [19-22], which is a low energy effective theory 
of QCD describing the dynamics of Nambu-Goldstone(NG) bosons, plays an important 
role to extrapolate physical observables such as NG boson masses and decay constants 
calculated in lattice QCD simulations at heavier quark masses to the physical quark mass 
point. Furthermore, not only the quark mass dependence but also the scaling violation 
for the Wilson quark action are described by the Wilson ChPT (WChPT), which includes 
effects of non-zero lattice spacing a [23-25]. Since 0(a) contribution can be absorbed into 
the quark mass term m, an inclusion of the lattice spacing effect in the ChPT is rather non- 
trivial, so that 0(a?) terms dominate in the small quark mass region. For QCD with the 
ordinary Wilson quark, several physical observables have been calculated in the WChPT 
at the next leading order (NLO), while only the leading order (LO) results exist in the 
WChPT for the tmlQCD 1 2 . In this paper, we therefore present results of pion masses 
and decay constants in the WChPT at the NLO, including 0(a 2 ,am) terms in the LO 
Lagrangian. 

In Sec. 2, we consider the power counting scheme in detail and give results at LO. In 
Sec. 3, we calculate the one loop contribution for pion masses and decay constants, and 
show that divergences of these quantities at one loop can be renormalized by the NLO 
counter terms. We present the pion masses and decay constants at NLO arbitrary value of 
the twist angle. In Sec. 4, we consider the automatic 0(a) improvement at the maximal 
twist, using the NLO calculation. We summarize our paper in Sec. 5. Details of NLO 
calculations are given in appendix. A. 

2 Analysis at "leading order" 
2.1 Power counting and Lagrangian 

In the ordinary continuum ChPT [19-22], the quark mass m, equivalently the meson mo- 
mentum p 2 , is considered as the expansion parameter, so that 0(M) terms consist of the 
LO Lagrangian where M = m or p 2 . At the NLO order, local counter terms of 0(M 2 ) 
cancel the divergences of one loop contributions generated by the LO terms. 

In addition to M, the lattice spacing a appears as the expansion parameter in the chiral 
perturbation theory for lattice QCD. Due to the explicit breaking of the chiral symmetry 
in the Wilson quark action, the lattice spacing effects start at 0(a). Therefore it is natural 
to treat 0(a) term as the LO contribution such that a ~ M. In the WChPT, this LO 0(a) 
term can be absorbed into the mass term by 

m — > fh = m + 0(a), (2-1) 

lr There exists the NLO ChPT calculations without 0(a 2 , am) terms for Nf = 2 [26-28], N f =3 [26, 28], 
and N f = 2 + 1 + 1 [30]. 

2 In ref. [31], NLO chiral log terms have been calculated with a similar power counting scheme, but in 
this paper we calculate and consider the renormalization at NLO in detail and the general case. 
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tree level 


LO 


0{M) 


SLO 


0{a 2 ,aM) 


one loop 


NLO 


0(M 2 ) 


NSLO 


0{a 2 M,aM 2 ) 



Table 1. Power counting scheme in this paper. We treat O(M) (LO) and 0(a 2 ,aM) (SLO) as 
the tree level Lagrangian. We introduce NLO and NSLO terms as local counter terms to cancel 
divergence of one loop contributions. 

so that no extra contributions to the continuum ChPT appear in the WChPT at this order. 
This power counting, however, is not correct in the region where rh is small, since the NG 
boson mass m n at the LO is given by 

ml = 2Bm, (2.2) 

which becomes tachyon for negative rh 3 . One has to add 0(a 2 ) terms to the LO Lagrangian 
for the stability of the vacuum to avoid the appearance of the tachyon [23, 25]. This LO 
Lagrangian is still insufficient due to the following reason. While 1-loop contributions from 
this LO Lagrangian generate 0(M 2 , Ma 2 ,a ) terms where M = rh or p 2 , terms odd in a 
such as O(Ma) are never generated from 1-loop contributions. Hereafter we ignore 0{a A ) 
terms since they are small in present lattice QCD simulations. 

From the above consideration, a physical observable X in general depends on both 
quark mass rh and a as 

X(rh,a) = rhX (m) + rhaX 1 (rh) + a 2 X 2 (rh) + O(a 3 ), (2.3) 

where we ignore small a 3 or higher contributions. We therefore set up the WChPT to 
calculate Xo(m), Xi{rh) and X2(rh) order by order in rh. For this purpose, in addition 
to leading order O(M) terms, we consider 0(a 2 , Ma) terms as the tree level Lagrangian. 
We call 0(a 2 , Ma) terms the sub-leading order (SLO) terms to distinguish them from 
LO terms. The one-loop contributions generate an additional factor of M to the tree 
level Lagrangian, LO plus SLO. Divergences from these one-loop contributions must be 
canceled by NLO 0(M 2 ) and NSLO 0(a 2 M,aM 2 ) terms. We stress again that higher 
order contributions such as a 2 M 2 are neglected here. 

For the twisted mass fermion [1, 2], there are two mass parameters, untwisted quark 
mass rh and twisted quark mass pi, which are denoted as m = rh, [i. Therefore M in this 
case represents m, /i or p 2 . In table 1, we summarize our power counting scheme. 

3 In the continuum ChPT, no tachyon appears for all m, since rh — ► \m\ in the mass formula, thanks to 
the chiral symmetry, which is absent in the WChPT. 
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2.2 Tree- level Lagrangian 



The tree-level effective Lagrangian for Nf = 2 tmlQCD, which include LO and SLO, is 
given by 



Cslo = W45<D A4 E^St)((S - Sq)^ + A(E - So) 1 ") 

- w 68 <s x t + x st)(sit + 1st) - w^<sit + i s t) 2 , 



(2.4) 



(2.5) 



where -D^S is covariant derivative with the left and right source current 1^ and r M defined 
by 



DfjTi = d^T, - il^T, + iEr^, 
= 8^+1%%-™^. 



(2.6) 
(2.7) 



After the construction of the Lagrangian, spurion fields x an d A should be set to 

X -»• 2S M = 25 (m + ir 3 /i), ^ 2W a. (2.8) 

The coefficients /o, Z>o , ^0,45,68 and Wg 8 in eqs. (2.4), (2.5) and (2.8) are the low energy 
constants, and their dimension are [/o] = [Bo] = 1, [Wo] = 3, and [W^s] = [W&s] = 0. 
The SU(2) matrix fields E(x) for NG bosons is defined by 



— ^o ^ph^O ' 



E ph = exp 



1 fo 



(2.9) 
(2.10) 



where So is a vacuum expectation value of E(x), r a is the Pauli matrices and 7T a (x) is the 
pseudo scalar NG field. 

As already mentioned, the O(a) term, (Y>A^ + AT^) can be absorbed to 0(m) term, 
(Ey^ + yE 1 ") by the replacement that 



2B m = 2B m + 2W a. 
Replacing the mass parameter m with shifted mass rh, we obtain 



fo 



£lo + slo = f [l + ^((E - So) + (S - Eo)^] (Z^EZ^St) 



/o 2 



^(E + Et)-^(E-EV) 



+ ^oggom - c 2 a (£ + E+)2 _ c^z^ + ^ _ 



16 



16 



where the coefficients c's are defined by 



C = 32^45 ^j, 



W 2 

C2 = -64(^8 -W 68 )-f, 
/o 



c 2 = 32W 6 



Wo 

• f2 ' 
JO 



(2.11) 



(2.12) 



(2.13) 
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2.3 Gap equation and pseudo scalar meson mass 

We first determine the vacuum in the Lagrangian (2.12). Parametrizing the vacuum ex- 
pectation value of £ as 

So = ex.p[i(pT 3 ] = cos 4> + it 3 sin </>, (2-14) 

the vacuum energy becomes 

V(4>) = — /q [2i?ofricos 4> + 2Bq{j, sin (j) — (c2<2 2 — C2(i2B§m) cos 2 cf) + c 2 a2i?oAi cos 4> sin 4>]. 

(2.15) 

The vacuum expectation value corresponds to the minimal point of this vacuum energy 
(2.15) (fro, which is determined by solving the gap equation, 

2Bqiti sin <f>Q — (c2d 2 — C2d2Bom) sin 2<f>o = 2B§\x cos 4>o + C2a2Bofx cos 2<f>Q. (2.16) 
We next expand the Lagrangian in terms of the component fields ir as 

£LO+SLO = ^LO+SLO + ^LO+SLO H ' ( 2 - 17 ) 

where 'Clo+slo represents the 0(iT n ) terms in the tree-level Lagrangian. With the com- 
ponents fields 7r a , the 0(it 2 ) Lagrangian ^q + slo * s Si ye by 

^cUslo^W + ^K) 2 - (2.18) 
where the tree-level pion mass m 2 a is written as 



[ (ml) 2 = m 2 + Am 2 (a = 3) 

to 2 = 2B$m! — 2c2d 2 cos 2 0o + 2c.2a(2B§m!) cos 0o, (2.20) 

Am 2 = 2c 2 a 2 sin 2 (f> + 2c 2 a(2B n') sin O - (2.21) 

Here m w = m^ (ml) denotes the charged (neutral) pion mass, and the short-handed 
notation for mass parameter is given by 



m \ I cos 0o sin <po \ I m 
H* I - sin (f) cos (f> \ fi 



(2.22) 



Using this notation, the gap equation (2.16) is written as 

2Bofi = — C2d 2 sin 20 + C2a(2Bom' sin0 — 2Bq\J cos 4>). (2.23) 

Using this recursively, ay! terms is found to be be of higher order, so that we can replace 
the /i' term with a 2 , am' terms. 

We now discuss the lattice cutoff a as well as quark masses m, /i dependences of vacuum 
condensation cos </>o and pion masses. Let us consider the case that the coefficient C2 = 0. In 
this case, there are two possible phase diagrams[24], depending on the sign of a coefficient 
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— 2 

2B m/c 2 a 

Figure 1. Vacuum angle cos </>o as a function of m for C2 > and 2Bo{i/c2a 2 — 0, 1, 2, 3. 




-6 -4 -2 2 4 6 -6 -4 -2 2 4 6 

2B Q m/ c 2 a 2 2B m/c 2 a 2 



(a) Charged pion mass (m^) 2 . (b) Neutral pion mass (m^) 2 . 

Figure 2. Pion masses as a function of m for C2 > and 2Bo/j,/c2a 2 = 0, 1, 2, 3. 

C2- Solving the gap equation (2.16) and minimizing the potential energy (2.15), we obtain 
the dependence of cos </>o on the untwisted quark mass rh. For C2 > 0, figure 1 shows the 
form of cos^o as a function of 2BQifi/c2a 2 at 2Bq^/c2Cl 2 = 0, 1,2, and 3. Corresponding 
pion masses are shown in figure 2. In the untwisted theory (// = 0), there are second order 
transitions at 2Borh/c2a 2 = ±2, as shown by the kinks in figure. 1 and by the vanishing 
pion masses in figure 2. The parity-flavor breaking phase, defined by the condition that 
cos 0o ^ ±1, lies between these two second order phase transition points, and two Nambu- 
Goldstone bosons associated with the flavor breaking appear in this phase [32-34]. Once 
2Bofi becomes non-zero, however, the transition turns into a crossover, and pion masses 
always stay non-zero due to the explicit flavor breaking by non-zero /x. For C2 > 0, the 
charged pion is heavier than the neutral pion, due to the 0(a 2 ) effect. If we change the 
value of the twisted mass continuously at fixed \2Borh/c2a 2 \ < 2, there appears the first 
order phase transition while crossing the parity-breaking phase at [i = 0. 

For C2 < 0, figure 3 shows cos^o and pion masses as a function of 2Borh/c2a 2 at the 
several fixed values of 2Bq[i. Figures 3-a and 3-b show the results for untwisted theory at 
fi = 0. The condensate cos^o jumps from Eo = 1 for 2Bqiti > to So = —1 for 2B$m < 0. 
This is the first order transition without the flavor breaking, so that all pions remain 
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(a) Vacuum angle at 2Bq^ = 




(b) Pion masses at 2Bq^ = 



2B m/c 2 a 




2B m/c 2 a 



(c) Vacuum angle at 2Bqh = C2<2 2 




(e) Vacuum angle at 2Bq[i = 2c2a A 




(g) Vacuum angle at 2Bqh = 3c2« 2 



(d) Pion masses at 2B§[i = 02a 1 
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(f) Pion masses at 2Bq[i = 2c2a 2 
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2B m/c 2 a 2 

(h) Pion masses at 2Bq[i = 3c2a 2 



Figure 3. Vacuum angle cos^o and pion masses as a function of rh for C2 < at 2£>o/x/ c 2 a2 
0,1,2,3. 



massive and degenerate. The effect of non-zero twisted mass fi generates a non-zero value 
of t 3 component, sin <fto = yl — cos 2 4>. There still remains the first order phase transition 
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(a) Vacuum angle for C2 > 
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(e) Neutral Pion mass for C2 > 
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(d) Charged pion mass for C2 < 
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(f) Neutral Pion mass for C2 < 





c 2 a=0 




c 2 a=0.2 




c 2 a=0.4 ; 



Figure 4. C2<2 dependence at 2B^ji = 



at which sin</>o flips sign between ±(1 — \2Bq^l/2c2Cl 2 \)- The first order phase transition 
disappears at \2Bq^/2c2Cl 2 \ = 2 and it turns into a cross-over at \2Bq^jl/2c2Cl 2 \ > 2. Due to 
the explicit flavor breaking, the neutral pion(-7r 3 ) is lighter than charged pions(7r ± ). 

These qualitative features remain true even for the case that £2 7^ 0. In figure 4(Left) 
we show results for C2 > at /i = 0. In this case the 2nd order transition points of the 
parity-flavor breaking phase move toward more negative values as the c^a increases. Pion 
masses are degenerate and their slopes increase (decrease) in the right (left) of the phase 
transition points in the symmetric phase, while charge pions become massless NG bosons 
and the neutral pion is massive in the broken phase. On the other hand, in the case that 
C2 < and fi = 0, the first order phase transition line does not move at all and the pion 
mass is still degenerate, as shown in figure 4(Right). The slope of both pion masses increase 
(decrease) in the positive (negative) rh region. 
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2.4 Vertices 



Since the twisted mass term /i explicitly break the parity symmetry, terms including odd 
powers of pion fields, ^^o+slo w ith n = 1,2,---, exist in addition to those with even 
powers, ^lo+SLO- F° r 71-3 terms, one finds 



^(3) 

£ LO+SLO = 7T f Md^a) 



H — \2Bqh' + &C2a 2 cos 0o sin 0o — \c2a(2B§m! sin 0o — 2B§y! cos 0o)]7r 2 7r3, 

6/o 



(2.24) 



while for the 0(ir 4 ) Lagrangian we have 

,(4) _ 1 



^LO+SLO ~~ g p 
1 



Kad^a) 2 ~ n + 7~C()0 COS (f)Q^J ^{d^af 

■[2Bom' — 8c2<x 2 cos 2 0o + 8c2a(2Bom!) cos </>o](vr 



24/ 2 ' 

- -^72 [8c 2 a 2 sin 2 0o + 8^0(2^') sin o ]vr 2 (^ 3 ) 2 . (2.25) 
2.5 Axial current and "decay constants" in twist base 

The quark bilinear operators are defined through the derivatives of the Lagrangian with 
respect to source terms as 

S» = W=^A J- (2.26) 

V = ^T* = 5(^ + ^) A A S = ?7,7 5 ^ = |(^-^)£. (2.27) 
Applying these derivatives to the tree-level Lagrangian (2.5), we obtain the axial current 



as 



A^ = -{r a (Ed^-^d^)) (2.28) 

(2.29) 



i/od^a cos 0o (« = 1,2) 



[z/od/i^a (a = 3) 

where in second line, we expand the £ in terms of pion fields 7r a up to the 0(ir) order. 
Since the pseudo scalar "decay constant" fps in twist base is defined from the expectation 
value of A"; between the vacuum and one NG boson state as 

(0|^(x)K(p)} = f a PS P,e-^5 ab} (2.30) 

we obtain 

n s = (f"*° ;° = 1 ; 2> (2.3D 

l/o (a = 3) 
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at this order. Note that this decay constants is in twist base, and we need the vector 
current for calculate it in physical base. 

The 7r 3 terms in the axial current ^4", needed at the NLO calculation, is given by 



(-4 



1,2n(3) 
M JLO+SLO 



2 
3/f 



1 



cq a cos <pQ 2 



2/ 2 



-IT d pit a 



cos <po 



co asm (p Q 
ij 75 ^3(7^0^3 - vr 3 d M 7r a J 



f- 



2 2 1 2 , Co a cos i 



-7r d M 7r a 



(2.32) 
(2.33) 



3 "next leading order" analysis 

3.1 NLO and NSLO Lagrangian 

The NLO and NSLO Lagrangians are given by 



£nlo = -O 



-^NSLO £-'ap 2 m,am 2 ^a- 



p^,a^m- 



(3-1) 



We here do not include the 0(p 4 , ap ) terms, since they do not contribute to pion masses 
and decay constants at this order. The 0(p 2 m,m 2 ) terms are the ordinal NLO terms in 
continuum ChPT [19-22] and are given by 



£p 2 m,m 2 



(3.2) 



In 0(p 2 m,m 2 ) terms, we follow the notation of [28]. The relation of the low energy 
constants in the notations between [21] I and [28] L is given by L45 = I4/8, Lqs = (^3+^4) /16. 

The 0(ap 2 ,am 2 ) and 0(a 2 p 2 ,a 2 m) terms describe the lattice artifacts and they are 
constructed as usual using the spurion analysis as 



J ~'ap 2 m,am 2 



+ ^ 6 (i x t + xi t )(Sx t + xS t ) 

+ V 23 (D„EA^D^ + D^xD^A) + V3a(^ 2 Sit + AD 2 E+)(Ext + x£ f ) 
+ V 3m (D 2 ^ + X-D 2 S f )(Sit + ASt), (3.3) 



^a 2 p 2 ,a 2 m 



+ x 4 (sit + ist) 2 (sx f + + x 5 (iit) (Sx f + xs f ) 
+ x 6 (ix t + xi t )(si t + ist) 

+ ^(^Sit^Sit + D^AD^A) 

+ X 3 {D 2 T,A lf + iL» 2 S t )(Si t + is f ), (3.4) 
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where the Vi(i = 1, . . . , 6), Vi(i = 23, 3a, 3m), Xi(i = 1,...,6), and Xt, (i = 23,3) are 
additional low energy constants. 

Here we introduce the nineteen parameters in NLO Lagrangian but our aim is not to 
determine these parameters. We instead investigate quark and lattice spacing dependences 
of physical observables using these parameters. In the physical observable, these parameters 
always appear as the some linear combinations, so that the number of independent fit 
parameters are much smaller. 

Expanding £ in terms of component fields tt, 0(tt) and 0{ir 2 ) terms of NLO+NSLO 
Lagrangian ^nlo+nslo and ^nlo+nslo are S iven b y 



^NLO+NSLO - CnLO 71 ^, 



-NLO 
•'-NLO+NSLO ~~ 

The coefficients Cnlo's are given by 



NLO^M 



{d^a? + C'NLO^a • 



(3.5) 
(3.6) 



/( 



Cnlo = — -{2B m'){2B fi')4L 68 - a{2B m') 2 (8V^ + sin</> 
i L 

+ a(2B m')(2B fi')(16Vl + 2Vq) cos <p - a(2B fJ,') 2 (V^ + 2Vq) sin^ 
- a 2 (2B m')(8X A + X' e ) sin20 o 

+ a 2 (2B 0f i'){(AX' A + X' 5 + X' 6 ) + (AX' A - X' 6 ) cos20 o }l , (3.7) 



n 



2p,a 
NLO 



2B m'L 45 + a ^2B m'(4V{ + V 2 ' - V^/2) cos (p - 2B^{y! 1 - V^/2) sin O } 
+ a 2 1{2X[ + X' 2 /2 - X' 23 /2) + 2X[ cos2<fo} 



a(2B 0l i')(2Vi + ~ WL ~ ^L) Wo 
-a 2 (X' 3 + X' 23 /2-X' 3 )(l- cos20 o )}<5 a3 



(3.8) 



C 



2,a 
NLO 



-{2B m') 2 2L m + a{(2B m') 2 (12Vl + V^/2 + V£) cos<^ 



16 

I 2 
Jo L 

- (2B m')(2B fj,')Ve sin^o + (2B fJ,')V^/2cos « 
+ a 2 (2B m') {{QX A + X' 5 /2 + X'J2) + {6X' A + X' 6 /2) cos 20 o - a 2 {2B [i')X' e /2sm2(j) } 
+ {(2B fi') 2 2L 68 - a8y 4 '{(25 o /x , ) 2 cos0o + 2(2B m')(2B 0f i') sin <fo} 



4XV {(2J3 m / )(l - cos 20) - 2(2B Q n') sm2<f> }}5; 



3« 



(3.9) 



where we use the normalized coefficients V/ = 2W V5, V( = 2W Vi, X[ = (2W Q ) 2 Xi and 
X[ = (2Wo) 2 Xi. Note again that, since the twisted mass term explicitly breaks the parity 
symmetry, odd power terms appear. In particular, the 0(ir) terms, together with ldoop 
contributions, must be canceled by the redefinition of the vacuum angle 4>q as 



4>R 



+ A0, 



(3.10) 



where A(p is the NLO contribution. 
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3.2 The vacuum expectation value and pion masses 

Due to the presence of three-point vertices (2.24), the tadpole diagrams £^ contribute 
to vacuum expectation value at 1-loop, 



4loop = f( A ± / ± + AoJ °) 7r3 ' ( 3 - n ) 

Jo 

A± = C2a 2 sin 2<f)Q + (co — C2)a(2Bom') sin (po, (3-12) 
3 1 

^-o = 2 C2fl2 sin2( Ao + -^(cq - 3c2)a(2B m') sin0 o , (3.13) 

where I a is given by 

Ia = I (2vr)V + K0 2 ' (3 ' 14) 
See the detail of the calculation in Appendix A. This contribution diverges and therefore 
must be renormalized by O(vr) terms in ^nlo+nslO' *°S e ^ ner W1 ^ n the redefinition of the 
vacuum angle 4>r. Explicitly the renormalization condition becomes 

'(i) , A 1 ) 



1 



,qj>LO+SLO 



^ + £iLp+ALO + NSLO = ; (3-15) 

0=00 



which leads to the renormalized vacuum angle (j)R = 4>q + Acfi as 
1 



W) 2 



A±L± + A L - -2 {a(25 m') 2 Ci + a 2 (2B m f )C 2 } 
Jo 



(3.16) 



The first term is the tree-level contribution for vacuum angle, while the second and third 
terms are chiral logarithm contributions from the charged and neutral pion, respectively. 
The last two terms are the polynomial contributions from the NLO Lagrangian. The chiral 
logarithm for the pion loop is defined by 

io .«£ log (^y, (3. 17) 

and the coefficients of the chiral logarithm, A®, are given in eqs.(3.12) and (3.13). The coef- 
ficients of the NLO polynomial term, C\ and C2, are some combinations of the renormalized 
NLO low energy constants. 

In the similar way, we obtain the result for the pion masses as 

+ {™lB a b +a 2 Qt)L b 
b=±,0 

1 f\ 

- -2 {El (2B m') 2 + E^a(2B m') 2 + E%a 2 (2B m')} . (3.18) 
Jo 

This is one of the main results of this paper. The first term is the tree-level form for pion 
mass but with the vacuum angle obtained at NLO. The second line represents the chiral 
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logarithm while the polynomial terms for the NLO Lagrangian are given in the last line. 
Each coefficient of the chiral logarithms is given by 



B± = a(2c - 4c 2 ) cos (j> , 
B°± = l + a(2c -2c 2 )cos<fo, 



Bq = - + o(c - c 2 ) cos (f) , 



Bq 



2 + a ( c o ~ 3c 2 )cos</>o, 



Q± 



4c 2 cos 



Q°± = 2c 2 cos 2 



<3o = C 2 cos2 00 

<5o = C2(3cos 2 i 



4 sin 2 0o )• 



(3.19) 

(3.20) 

(3.21) 
(3.22) 



and coefficients of polynomial terms Ei(i = 1,2,3) are combinations of the renormalized 
NLO low energy constants. Our results have not only the same pion mass contributions 
for the chiral logarithms in [31], but also the additional 0(am) terms from SLO. Note that 
0(1) terms in eqs.(3.19) and (3.20) of course agree with continuum results [19, 20]. 

3.3 "Decay constants" 

Let us consider the NLO axial currents. Using eq. (2.27), we obtain the "NLO" (=NLO+NSLO) 
terms of the axial currents for a = 1,2, 



16* 



NLO 



0, 



. -'Ii'^a 
JO 



{(2B m')L 45 + a{2B m')(W{ + V 2 ' - V 2 ' 3 /2 - 2V£ m ) cos O 

- a(2B ^')(V 2 ' - Vlz/2 - 2Vi m ) sin0 o 

+ a 2 {2X[ + X' 2 /2 - X' 23 /2 - 2X' 3 + 2X[ cos 2cb )} COS (f>Q 



- a(2B m')2V 3 ' a 



(3.23) 



and for a = 3, 



16* 



(4)nlo = -T d ^z [{(2B m')L i5 + a(2B m')(W{ + V± - U 2 ' 3 /2 - 2V 3 ' a - 2F 3 ' m ) cos , 



- a{2B^')(Vi + 2V 3 ' + Viz/2) sin O 

+ a 2 {2X[ + X' 2 /2 + X' 3 -X' 3 + (2X[ - X' 3 - X' 23 /2 - X' 3 ) cos 2^ } 



(3.24) 



From these results, we calculate the "decay constants" for the pseudo-scalar meson at 
"NLO", which are given by 



fp 



PS, NLO 



PS,LO 



+ fpS,LO [ ^2 



^b 



6=±,0 



16 



+ -J {{2B m ! )Hl + a(2B m')H^ + a 2 H%) 
Jo 

+ Af ± (l-S 3a ). 



(3.25) 



This is also one of the main results of this paper. The first term is the tree-level contri- 
bution, which includes the vacuum renormalization effects at "NLO". The second term 
represents the chiral logarithms and the polynomial terms for the NLO Lagrangian are in 
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the second line. In addition, the "decay constant" for the charged pseudo scalar has the 
additive renormalization term in the last term. Coefficients of the chiral logarithms are 
given by 

F± = ~ (1 + c a cos 0o ) , Fo = ~ U + ^c a cos O 1 , (3-26) 

F± = - (l + ^c a cos , F ° = -^c a cos 0o, (3.27) 

and the formula for the coefficients of the renormalization terms Hi(i = 1,2,3) and Af± 
are given in Appendix A. Note again that 0(1) terms in eqs.(3.26) and (3.27) agree with 
continuum results [19, 20]. 

3.4 Consistency 

Divergences generated by 1-loop contributions must be removed by the "NLO" terms. For 
this purpose we define the renormalized NLO and NSLO low energy constants as 

Li = L[( MC hPT) + (3.28) 
V( = Vf(MchPT) + ^R, (3.29) 
Xi = Xj(f, Ch PT) + ^R, (3.30) 

where h, renormalize parameters, the argument /ichPT represents the renormal- 

ized scale and R = 0(l/e) is introduce to cancel 1-loop divergences in the dimensional 
regularization. We give the detail formulation in appendix. A 

The renormalization conditions for the vacuum expectation value leads to 

8(8w 4 + v 5 + 4/ 68 c 2 ) = 3c - 5c 2 , (3.31) 
8(8x4 + x 6 - 4/ 68 c 2 ) = 5c 2 . (3.32) 

In order to cancel the divergences in 1-loop diagrams for pion masses, coefficients 
should satisfy the following conditions. 

16(/ 45 - 2l m ) = 1, (3.33) 

16(12« 4 + v 5 /2 + v 6 + 4wi + v 2 - V23/2 + 2c 2 Z 45 ) = -6c 2 + 6c , (3.34) 

16(6x4 + x 5 /2 + x 6 /2 + 2xi + 2:2/2 - x 23 /2 - c 2 / 45 ) = 4c 2 , (3.35) 

16(6x4 + x 6 /2 + 2xi - c 2 Z 45 ) = 2c 2 , (3.36) 

16(4x4 - x 3 - x 23 /2 + x 3 - c 2 / 45 ) = -6c 2 . (3.37) 
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Similarly from the decay constant, we have 



4Z 45 = 1, (3.38) 

I6v 3a = c , (3.39) 

8(2vi + v 2 /2 - £ 23 /4 - 2v 3m ) = 2c 2 + 7c /4, (3.40) 

16xi = -3c 2 , (3.41) 

16(x 2 /4-x 2 3/4-2S 3 ) = 2c 2 , (3.42) 

l6(x 3 /2 + x 23 /4 + 3x 3 /2) = -c 2 . (3.43) 

A solutions to these conditions can be obtained as 

k 5 = \, hs = ^, (3.44) 

3 i 23 7 

xi = -t^c 2 , x 2 - x 2 3 = -c 2 , x 3 + — = — c 2 , (3.45) 

lb 2 lb 

x\ = gC 2 , x 5 = ^c 2 , x 6 = 0, (3.4b) 

x 3 = ~C2, (3.47) 
lb 

16«i + 4v 2 - 2£ 23 - lb£ 3m = 7 -c + 2c 2l (3.48) 

b4u 4 + 8v 5 = 3c -8c 2 , (3.49) 
b4v 4 + 8v 6 + 32£ 3m = -2c -7c 2 , (3.50) 

ha = ^. (3.51) 
Note that results of Z45 and ^8 agree with the continuum ChPT [19, 20]. 

4 maximal twist 

One of the advantage of the tmlQCD is the automatic 0(a) improvement at maximal twist 
[3, 5]. The ETMC employs the PCAC quark mass to determine the maximal twist in their 
simulation, 

mpcAc(/i) - 2E^MPl) " (a ~ 2) - 

The above condition must be satisfied by tuning the untwisted mass m at each twisted 
mass /j,[5]. The untwisted mass which realizes the maximal twist condition is called the 
critical untwisted quark mass. This definition for the maximal twist is called the PCAC 
definition. In order to reduce the numerical cost, however, m for the maximal twist is 
defined at /U m ; n , the minimal value of the twisted mass employed in the simulations, such 
that mpcAC (thain) = 0. This is called as the fixed PCAC definition. In figure 5, both 
definitions are schematically drawn for c 2 < 0. In this section, we investigate the difference 
between two maximal twist conditions, using our WChPT analysis. 
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Figure 5. The chiral limit for c-i < 0, where the first order phase transition line exists. The 
untwisted mass m corresponding to the maximal twist depends on twisted mass fj, in the PCAC 
definition, while m is constant in the fixed PCAC definition. 



4.1 PCAC quark mass 

The PCAC quark mass can be written in the WChPT as 

fpsmj 

mpcAC = (4.2) 

where fps, m -w denote the decay constant and mass for the charged pion, and we define 
Z(p) as 

(0\P a (x)\7T(p))=iZ(p)e-^. (4.3) 

At "NLO", we have 

/p5,nlo = /ps.loI^-^+a^ (! + A /) + A/±, (4.4) 

ml tNh0 = 7< LO U^ o+A0 (! + A -)> ( 4 - 5 ) 

(1 + A Z ), (4.6) 

where the Aj(i = f,m,Z) denote the NLO contribution corresponding to fps-, m -K and 
Z(p). The PCAC quark mass at the NLO is given by 



mpCACNLO = m PCAC,Lo|0 o _>0 o+A (1 + A/ + A m - A Z ) , 

A/± (4.7) 



9 

mz 



V,LO 



2Z LO (0) 
where 



o -></>o+A0 



mpcACLol^^ = 2Bo(1 + g 2< £ 8 (^ o + A0)) cos ^° + A ^ ^ 
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4.2 PCAC definition 

The maximal twist condition in the PCAC definition leads to 

cos 0o = 0{am') = 0{a^). (4.9) 

This condition simplify the gap equation (2.16), and the critical untwisted quark mass as 
a function of the twisted mass [i is given by 

2B m = -c 2 a2B fi + 0{afi 2 ). (4.10) 

Eq. (4.10) shows that the critical untwisted quark mass rh depends linearly on the twisted 
quark mass /i. Using eqs. (4.9) and (4.10), the charged pion mass at the NLO is given by 



(^)nlo = ( m t)lv 

L Jo 



1 + ^o-^(«)!o^ + « 2 ^) 



(4.11) 



where, the NLO polynomial coefficients L r c and X r c are given by 



L r c = Ll 5 -2L r 68 , (4.12) 

x r c = (x r 5 /2 + xy2-xy2), (4.13) 

the chiral logarithm Lq is defined by eq. (3.17), and the LO pion masses are simplified as 

(m±) 2 = 2B », (4.14) 
(m°) 2 = 2B 0f i + a 2 {2c 2 + 2c 2 ,2£ /i). (4.15) 

Eq. (4.11) shows that the charged pion mass is 0(a) improved: the lattice spacing cor- 
rections start at 0(a 2 ). It is important to note that the charged pion mass at the NLO 
contains the chiral logarithm only from the neutral pion loops but not from the charged 
pion loops at the maximal twist, as in the continuum ChPT. 

Note that the maximal twist by the PCAC quark mass and the charged axial currents 
give the same constraint for physical observable in this WChPT analysis. 

4.3 Fixed PCAC definition 

For the fixed PCAC definition, the maximal twist condition and the critical untwisted 
quark mass are given by 

cos 0o = 0(a/i min ) (4.16) 

2B m = -c 2 a2B ii Tain + 0(ap 2 min ), (4.17) 

where /U m ; n denotes the minimal value of the twisted mass, used for the fixed PCAC defi- 
nition. Eq. (4.17) shows that the critical untwisted quark mass is constant and does not 
depend on the twisted mass /z. Using eqs. (4.16) and (4.17), the charged pion mass at the 
NLO is given by 



(^)nlo = (™t)l 



0o->0o+A0 



l + lL -^((mt)l o L r c + a 2 X r c ) 
z Jo 



(4.18) 
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where the coefficients L r c and X r c are defined in eqs. (4.12) and (4.13) and the pion masses 
are given by 




</>o-*£o+A</> 



2B n, 

2B n + a 2 (2c 2 + 2c 2 2 2B n min ). 



(4.19) 
(4.20) 



The charged pion mass (4.18) for the fixed PCAC definition has the same functional form 
as (4.11) for the PCAC definition. The fixed critical untwisted quark mass effect, however, 
appears in the 0(a 2 ) terms for the LO neutral pion mass. 

It turns out that the charged pion mass is 0(a) improved for both PCAC and fixed 
PCAC definitions, and that the difference between these two definitions appears as the 
0(a 2 ) effects. In addition, contrary to the case at the non-maximal twist where the chiral 
logarithm of the charge pion mass contains effects from both neutral and charged pion 
loops due to the lattice artifact, it contains only the effect from the neutral pion loop for 
both definitions as in the continuum ChPT. 

5 Conclusions 

In this paper, we construct the Wilson chiral perturbation theory for the Nf = 2 twisted 
mass lattice QCD at the small quark mass regime such that m q ~ a 2 A 3 . In order to consider 
such a regime, we include 0(a 2 ) and 0(am) terms at the tree level as the sub-leading order 
Lagrangian, which induce the non-trivial phase structure and pion mass splitting at the 
tree-level. Using this effective theory, we investigate the pion mass and decay constant 
as a function of not only the twisted quark mass but also the lattice cutoff at the next 
leading order. Our main results are given in Eqs. (3.18) and (3.25). We also confirm that 
divergences from 1-loop contributions can be consistently removed by the next leading 
order Lagrangian. 

For the comparison of our results with data obtained by numerical simulations, we 
derive the twisted quark mass dependence of the charged pion mass at the maximal twist. 
As the definition of the maximal twist, we adopt two different definitions, the PCAS and 
the fixed PCAC definition, the latter of which is actually employed in the simulations. We 
have found that the charged pion mass is 0(a) improved, so that lattice spacing corrections 
start at 0(a 2 ) for both two definitions, and that the difference between the two definitions 
appears as 0(a 2 ) effects. In addition, it should be noted that the chiral logarithm in the 
charged pion mass comes from the neutral pion loop only, as in the continuum ChPT. 
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A detail of calculations 



A.l NLO terms 

From the spurion analysis, 0(ap 2 , am 2 ) terms are given by 

0(ap 2 m) (Sit + ist)(Sx f + x^)(D^D^), (Atf + x^)(D^D^), 

(D^A* + AD^)(D^ X ] + XD^), (D^AW^ + D^xD^A), 

V v ' 

Si 

(D 2 Eit + AD 2 Yt){Y,x} + xS f ), (Z^Ex 1 + x£» 2 St)(sit + Atf), (A.l) 



s 2 S3 

0(am 2 ) (Sit + ist}(Sx f + X^f, (xx*)^ + (i* f + x^Ptf + X^}- 

(A.2) 

Here Si, S2, and 53 terms come from the non-commutativity between the covariant deriva- 
tive Du and mass term x- Therefore they vanish for the degenerated Nf = 2 (untwisted) 
Wilson fermion case. In the L = r« = limit, these terms can be expressed by other 
terms, 

51 i ((^ t Sit + i^St}(^Sxt + x^St) - (i^t + x^X^E^St)) , (A.3) 

5 2 -(^Sit + i^St)(^Sx f + x^St), (A.4) 

5 3 -> -(^Sit + ia /x Et)(^Sx t + X^St). (A.5) 
In a similar way to 0(ap 2 m,am 2 ) terms, 0(a 2 p 2 ,a 2 m) terms are given by 

0(ap 2 m) (Sit + i S t) 2 (^S^St), (AA^)(D^D^), (D^A + AD^f, 

(D^Sit^sit + DutfAD^A), (L> 2 Sit + iD 2 st)(sit + ist), (A.6) 

S v ' V v ' 

S4 S5 

0(am 2 ) (Sit + i S t) 2 (Sxt + ftf), (AA^ptf + X^), (Atf + x#)pA* + ASfi), 

(A.7) 

where S4, S5 terms vanish in the untwisted theory. In the L = = limit, they become 
S 4 -> ^(^Sit + i^St)(^Sit + i^St) - (iit)(^Sa M St), (A.8) 
S 5 ^-(^Sit + ia M St) 2 . (A.9) 

A.2 Renormalization for a vacuum expectation value 

Since there exist three point functions in the Lagrangian (2.24), we have one loop diagrams 
for a vacuum expectation value. We obtain 



o 




a2Bom'(2co — 2S2) sin <pQ + 2c2a 2 sin 2(/)q 
Wo 

a2Bom'(co — 3c2) sin cpo + 3c2a 2 sin 2<pQ 



I±, (A. 10) 



2/o 



J , (A. 11) 
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where the single(double) line represents the charged (neutral) pion, I±/q denotes a contri- 
bution from a charged or neutral pion loop, which is given by 

d 4 p 1 



L 



±/o 



2 

±/0\2 



(27T;- p A + ( m 



16vr 2 



2 [m 
— + 7js - 1 + log 
e 4-7T 



(A.12) 



where e is the dimensional regulator and 'Je is the Euler-Mascheroni constant. One-loop 
effects to the vacuum expectation value are summarized as 



C 



(i) 

Hoop 



-L(A±I± + A I )7r 3 , 
Jo 



A± = c 2 a 2 sin 20o + (co — c 2 )a(2Bom') sin 
A 



3 1 

— c 2 a 2 sin 20o + — (cq — 2>c 2 )a(2BQm!) sine 



(A.13) 
(A.14) 
(A.15) 



Using the gap equation (2.23), on the other hand, contributions from £nlo+nslo & r e 
given by 



C 



(i) 

NLO+NSLO 



-— [a(2B n 
Jo 



f\2 



(8T4 + V5 + 4c 2 L 6 8)sin< 



+ a 2 (2B m'){8X i + A 6 - 4c 2 L 68 ) sin20o]7r 3 . 



(A.16) 



In order to cancel divergences from the 1-loop integral, low energy constants in £nlo+nslo 
need to be renormalized, and renormalized low energy constants are given by 

U ~ (A.17) 

(A.18) 
(A.19) 

(A.20) 



L i — ^(MCIiPt) + ^3^' 
V( = ^(MChPT) + 



32vT 2 



X i — -XK/^ChPT?) + 



32vr 2 



R, 



where the argument /UchPT is a renormalization scale and R is defined by 

2 



R 



log(4vr) +7^-1, 



which cancels the 1-loop divergence in dimensional regularization. To renormalize the 
vacuum expectation value correctly, we need renormalization conditions, eqs. (3.31) and 
(3.32) , and the renormalized Lagrangian for tt is written as 

l±<p -f Z- lloop i- -^-NLO+NSLO 



'(1) 



-C 



LO+SLO 



(m°) 2 / o A0 + fo(A±L ± + A L ) - — {a(2B m') 2 d + a 2 (2B m')C 2 } 

To 



where the low energy coefficients C\, C 2 are given as, 

d = {8V[ + VI + Ac 2 L r 68 ) sin 0o, 
C 2 = (8X1 + XI - Ac 2 L r 68 ) sin 20 o . 



(A.21) 

(A.22) 
(A.23) 
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A. 3 Renormalization for mass 



The 1-loop diagrams, which contribute to the pion mass term, can be classified into two 
types. The first type are made from a four-point vertex, while the second ones are made 
from two three-point vertices. As the first type, we have following diagrams. 



In the second type, there are three diagrams. 




In the same way as tadpole diagrams, the contributions of the first type diagrams are given 
by 



1 

2/? 



, [Zt{d^±f + (AM±) 2 4] J±, (A.24) 



-L [Z±{d^ ± f + (AM±)\l] J 0l (A.25) 



-L [Z ± (d,7r f + (AM° ± ) 2 4] J ± , (A.26) 



2 [Zl(d^f + (AM °) 2 7r 2 ] J . (A.27) 



where each coefficient is given by 



Z± = — ^ — c a cos 00) (AM±) 2 = m 2 j-^ + a(c - 4:C 2 ) cos0 o j + 4c 2 a 2 cos 2 O , 

(A.28) 

± 1 c o acos0 o /A , /f ±N2 2 J 1 , a(c - 2c 2 ) cos O \ 2 2, / A on\ 

^o = - 3 2 ' ' ' = m?r I 6 2 I C2 ° COS ^ ^ 

Z± = ~- c a cos 0o, (AM^) 2 = m 2 | i + a(c - 2c 2 ) cos0 o 1 + c 2 a 2 (2 cos 2 00 - ^ sin 2 o ) 



(A.30) 

c a cos 0o / A /if(h2 2 / 1 , a(c -6c 2 )cos0 o ) 2 2 

Z = , (AM ) =m n -[-- + } + c 2 a (3cos O - 4 sin <p 



(A.31) 



- 21 - 



On the other hand, the second type diagrams give SLO 2 effects as 




(SLO) 5 



SLO 2 . 



p2 _|_ m 2 p2 _|_ m 2 



(A.32) 



Therefore the second type diagrams are higher order than NLO and we do not consider 
them in this paper. We then obtain the renormalized Lagrangian for ir 2 at 1-loop as 



C 



(2) _ r {2) r (2) r (2) 

~ ^LO+SLO ^lloop ^NLO+NSLO' 

2 I /o 2 



C 2 ^o} (W + \ {(<) 2 + ^ + C 2 - } vr 2 . (A.33) 



From this result, the renormalized mass is given by 



a\2 



NLO 



:o 2 + 



,a\2 



,m~)- < 1 

a\2 



f 2 
JO 



+ Ci 



2,a 
NLO 



io 



2p,a 1 



'NLO J 

(, m 7r)Lo|0 o ^^ o+ A</, 

+ ]T (m 2 ^ + a 2 Q^)L b 



JO 

Mgl b 



+ 



ft 



+ C, 



NLO 

2,a 
NLO 



6=±,0 

16 



{Ef(2B m') 2 + E%a(2B m') 2 + E^a 2 (2B m')} , (A.34) 



o 



where coefficients of the chiral log terms are given in eqs. (3.19)^(3.22), while 
Ei = L\ h — 2Lg 8 , 

E% = (12Vf + V£/2 + Vg + 4V{ + V{ - V{ 3 /2 + 2c 2 L r 45 ) cos cf> , 
El = (6X1 + Xl/2 + Xl/2 + 2X1 + X r 2 /2 - X r 23 /2 - c 2 L r 45 ) 
+ cos 20 o (6X 4 r + X£/2 + 2X{ ~ c 2 ^4 5 ) 
- 5 a0 (4Xl - XI - X r 23 /2 + XI - c 2 Ll 5 )(l - cos 20 o ) 

for NLO low energy constants. 



(A.35) 
(A.36) 



(A.37) 



A. 4 Renormalization for decay constant 

Since we have 0(vr 3 ) terms (2.32, 2.33) in the axial current, the 1-loop contribution to the 
decay constant of the charged pion is given by 



<°KKW>lloop 



c 



f,a 

Hoop 



3/o 2 



+ 



Co a cos <po 



3/ 2 



2 \ (2 3coa cos i 
- + c a cos (p I L± + I - H - — 



COS 00 

L 



coasm <pQ 

I 2 
Jo 

COS (f>Q — CqClLq. 



(7I"37T 3 ) 



(A.38) 
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For the neutral pion, 



(0\A 3 Mp)) 



a 



Hoop 

/,3 
Hoop 



-foP»C f ' 3 



3/ 2 



+ 



Hoop' 

Co a cos 4>q 



•J Jo 



4 \ c a cos 4> T 
- + c a cos 0o I L± H L . 



(A.39) 



Therefore, we obtain the renormalized decay constant 

n -1/2 



(o|^2K(p)>| 



1 4- 4- r 2p 

jo 



pi a ^■"^NLO , 

(A.40) 



where the first term is a contribution for the renormalization of the pion field, and the 
second and third terms represent the 1-loop and the "NLO" contributions, respectively. In 
the same way as the vacuum expectation value and the pion mass, we obtain the decay 
constant at NLO as 



JPS,NLO — JPS,LO\ 

16 



+ fps,t l r b 

6=±,0 



\<t>=(j> Q +A4> 

+ -2 {(2B m')H? + a(2B m')H$ + a 2 H%} 
Jo 

+ A/±(l - 5 3a ). 

where coefficients for the chiral log terms are given in eqs. (3.26) and (3.27), while 
Hi = £45/2, 

H- = 2V{ + V{/2 - VJ3/4 - 2V£ m - 2V£J 3a , 
£T« = X{ + X 2 r /4 - A^'3/4 - 2X3 + X[ cos 20o 

+ (Xl/2 + X 2 r 3 /4 + 3A > 3 72)(1 - cos 2<P)5 3a 



(A.41) 



(A.42) 
(A.43) 

(A.44) 



for NLO low energy constants. Note here that an additive renormalization term for the 
decay constant of the charged pion exists: 



A/± = -f [-a(2B m')2V 3 r a + c aL .] 
We finally consider the PCAC mass, defined by 

EJl^(0\dai(x,t)\n( P ))^(7r(p)\P a iO)\0) 
2E,/7&(0|P a (^t)k(p))^(7r(p)|P«(0)|0) 



(A.45) 



"1PCAC 



fpsml 

2Z PS 



(A.46) 
(A.47) 
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Using the same spurion analysis as for the axial currents, the pseudo scalar density and its 
renormalization factor Zp S are given by 

P a = i/ o £o(l + c 2 acos0oK (a = 1,2), (A.48) 
Z£ s = f B (l + c 2 a cos O ) (A.49) 

at the tree-level. Therefore the LO PCAC quark mass becomes 
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